UNCLASSIFIED 


AD  NUMBER 

AD078606 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  Oct  1955. 
Other  requests  shall  be  referred  to  U.S. 
Army  Ballistic  Research  Laboratories, 
Aberdeen  Proving  Ground,  MD . 

AUTHORITY 

USARL  ltr,  24  Jan  2003 


THIS  PAGE  IS  UNCLASSIFIED 


p*^ 


mm 


« &  i 


clinical  Information 


Reproduced  by 

CUMENT  SERVICE  CENTER 

OU^OTT  BUILDING,  DAYTON,  2,  OHIO 


This  document  is  the  property  of  the  United  States 
'loveritment.  It  is  furnished  for  the  duration  of  the  contract  and 
shall  be  returned  when  no  longer  required,  or  upon 
recall  by  ASTIA  to  the  following  address: 

Armed  Services  Technical  Information  Agency,  Document  Service  Center, 

Knott  Building,  Dayton  2,  Ohio. 


rv, 

Nv 


GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
OH  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFINITELY  RELATED 


PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
EE3 K» -ISIBIIUTY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THK  FACT  THAT  THE 
MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
»  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
-LICATi  W  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 


ISON  Or  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
;  OR  SKLL  AliY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


REPORT  No.  961 


An  Analytical  Treatment 
Of  The  Problem  Of  Triangulation 
By  Stereophotograimmetry 


HELLMUT  H.  SCHMID 


DEPARTMENT  OF  THE  ARMY  PROJECT  No.  5B0306011 
ORDNANCE  RESEARCH  AND  DEVELOPMENT  PROJECT  No.  TB3-OS38  gC'wd 

DAI  !  ICTI  r  DCCCADfU  I  ADHPATAPICC  b 


BALLISTIC  RESEARCH  LABORATORIES 


ABERDEEN  PROVING  GROUND,  MARYLAND 


B  A  L  t.  I  STIC  'RESEARCH  LABORATORIES 

REPORT  NO,  961  " 

(<  -  '  .  ' 

Vi  /'  .1  "  ' 

October  1955 

"...  // 

.  ..  „  /• 


AM  AXilr^CAL  TREATMENT  of  tha  PROBLEM  of  TEUNQUUTION 
by  STERSOPHOTOORAlMRar 

Ha  limit  H.  Schmid 


Dapartmant  of  tha  Army  Pro j act  Mo,  5 BO 306011 
Ordnarca  Maaaaroh  and  Daralopmarrt  Projaot  Mo.  TB3-0538 


ABERDEEN  PROVING  HOUND, 


MARYLAND 


■■  ,  -  .  M  '■  ..  '  -  -  ?■  ,.  .. 

*  ’  v  > 

.. 

II 

.  '  '  '  -  *;?  : 

*!> 

„  ..  ...  "  ”  >t 

o 

if 

"  ,  / 

.  TABLE  of  doNTENTS 

/*  ‘  . 

-  \ 

> 

vs 

v  ■  ■  :  ■  1  ■  , 

- 

<•  /i)  ^  1 

Abstract  * 

vs. 

3 

■  I. 

«•  "...  ..  ..  ii' 

Int  induction 

$  . 

,v  II. 

Qeneral  remarko 

5 

" 

- 

III 

.  The  ^«jpm®trical  analysis  » 

6 

IV. 

The  least  equanss  adjustment 

18 

18 

A.  The  process  of  orientation 

1.  The  least  squared  solution  -  «  v. 

\\ 

no 

■  >*w  > 

2*  The  derivation  oi‘  the  coefficients  of  the 

29 

observational  equations 

♦  ; 

1. 

n 

B.  The  prooess  of  triangclation 

36 

1.  The  triangulation  as  part  of  the  prooess 

36 

" 

♦ 

of  orientation 

, 

2.  The  triangulation  as  an  independent 

37 

l 

i 

computational  procedure 

C.  The  dete  rod  nation  of  the  mean  errors  of  the  *- 

38 

i  . 

// 

...  v  t  *  *  ‘  . 

observations,  of  the  elements  of  orientation, 
and  of  the  triangulation  results 

i  \ 

(0) 

..  ( 

V. 

The  application  of  the  analytical  method  to  the  problem 
of  control  extension 

Uo 

VI. 

Concluding  remarks 

5U 

•» 

References 

% 

! 

H 

■  n  .. 

2 

-  '  ' ;  ^ 

‘\ 

'.  1 

BALLISTIC  RESEARCH  LABORATORIES  REPORT  NO.  961 


ESchmid/mal 

Aberdeen  Proving  Ground,  Md. 
October  1955 


AIT  ANALYTICAL  TREATMENT  of  the  PROBLEM  Of  TRI ANGULATION 
by  STERECPHOTOGRAMMETRY 


ABSTRACT 


An  analytical  solution  for  the  ba3ic  problem  of  triangulation  by  stereo- 
photogrammetry  is  derived.  The  most  general  case  Is  defined  as  -the  problem  of 
determining  simultaneously  the  orientations  of  the  two  cameras,  whereby  all  18 
unknown  parameters  of  the  orientation  are  considered  with  no  limitations  on  the 
camera  orientations.  The  process  of  triangulation  is  treated  as  the  computation 
of  spatial  coordinates  as  functions  of  the  elements  of  orientation  and  the 
corresponding  plate  measurements .  The  least  squares  solution  derived  is  based 
on  rigorous  mathematical  expressions  which  connect  the  plate  measurements  with 
the  unknown  parameters .  In  contrast  to  the  conventional  approach,  the  separation 
of  the  orientation  problem  into  the  two  phases  of  relative  and  absolute  orienta¬ 
tion  is  avoided.  This  analytical  solution  can  be  based  on  a  few  basic  theorems 
of  solid  analytical  geometry.  Because  the  observation  method  -  monocular  or 
stereoscopic  -  does  not  influence  the  formulas  expressing  the  rigorous  geometry, 
it  is  possible  to  nwke  use  of  absolute  control  points  which  axe  not  common  to 
the  area  covered  by  the  two  photographs  under  consideration.  Thus  more  favorable 
geometry  is  introduced  into  the  problem  of  the  double -point  intersection  in  space 

Ihe  least  squares  solution  derived  is  suitable  for  any  number  and  any  com¬ 
bination  of  absolute,  partially  absolute  ard  relative  control  points.  In  addition 
any  one  of  the  elements  of  orientation  -  including  the  base  line  components  -  may 
be  enforced  in  the  solution.  By  applying  the  concept  of  pseudo-residuals  and  by 
introducing  cross -weights,  it  is  possible  to  treat  the  least  squares  solution 
like  a  problem  involving  independent  indirect  measurements.  Furthermore,  iw  ic 
shown  that  the  normal  equation  system  can  be  formed  step  by  step.  This  method  has 
msrit  -hen  electronic  computers  are  used  since  the  number  of  points  carried  in 
the  solution  has  only  a  slight  effect  on  the  amount  of  memory  space  needed. 

The  introduction  of  rotational  auxiliaries  which  are  essentially  direction 
cosines  and  the  combination  of  these  with  the  plate  coordinates  as  linear 
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auxiliaries  render  the  coefficients  of  the  observational  equations  c"  partial 
differential  quotients  in  terms  attractive  for  an  analytical  treatn^nt.  The 
process  of  triangulation  is  treated  as  a  part  of  the  process  of  orientation  as 
well  as  an  independent  computational  procedure. 

A  special  chapter  deals  with  the  determination  of  the  mean  errors  of  the 
observations ,  of  the  elements  of  orientation,  and  of  the  triangulation  results. 

Finally,  the  application  of  the  proposed  analytical  method  to  the  problem 
of  control  extension  is  discussed  in  principle.  It  is  shown  that  the  mel hod  used 
on  the  universal  plotters  for  a  strip  triangulation  procedure  is  an  approximate 
solution  only,  because  it  is  based  on  incomplete  conditional  equations. 

The  rigorous  geometry  for  the  problem  of  extension  is  interpreted  ee  the 
condition  that  rays  originating  from  three  consecutive  camera  stations  have  to 
intersect  for  at  least  one  point  located  in  the  area  common  to  the  three  photo¬ 
graphs  under  consideration.  The  corresponding  conditional  equations  are  derived 
and  the  coefficients  of  the  corresponding  observational  equations  are  given.  It 
is  shown  that  it  is  now  possible  to  include  in  the  extension, models  which  are 
formed  by  the  combination  of  photographs  taken  at  every  other  camera  posit ion. 
The  thus  extended  base  line  provides  for  a  favorable  base -height  ratio  otherwise 
obtainable  only  by  convergent  photography. 

As  stated  before,  the  final  normal  equation  system  can-be  formed  step  by 
step.  Attention  is  called  to  the  fact  that  the  matrix  of  the  unknown  parameters 
is  filled  in  the  neighborhood  of  the  diagonal  only,  thus  making  it  possible  to 
use  an  iterative  subroutine. 

The  method  presented  for  a  strip  triangulation  is  useful  in  a  least  squares 
treatment  of  a  block  triangulation  also. 
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I.  INTRODUCTION 


Fhotcgrammetry  may  be  defined  as  the  science  in  which  geometrical  properties 
of  objects  are  analyzed  in  quantitative  terms  from  their  images  recorded  on 
photographs.  Stereophotogremmetry,  in  particular,  deals  with  the  triangulation 
of  pencils  of  rays  originating  from  two  camera  Btations  by  applying  the  technique 
of  stereoscopic  observation. 

The  process  of  triangulating  two  corresponding  pencils  of  rays  consists  of 
the  restitution  of  the  orientation  of  the  two  photographs  under  consideration, 
and  the  reconstruction  of  the  model.  The  analytical  equivalent  of  the  restitution 
of  the  orientation  is  the  determination  of  two  sets  of  nine  degrees  of  freedom, 
namely,  six  translations  and  three  rotations  for  each  of  the  two  camera  stations. 
Expressed  analytically,  the  reconstruction  of  the  model  is  the  process  of  tri¬ 
angulating  individual  points  whereby  each  triangulation  is  based  on  the  18  ele¬ 
ments  of  orientation  and  four  corresponding  plate  measurements. 

Approaching  the  problem  analytically,  for  the  general  case  the  18  elements 
of  orientation  are  obtained  from  a  least  squares  adjustment  based  on  a  sufficient 
number  of  conditions  of  intersection  (relative  orientation)  and  on  independently 
established  control  points  (absolute  orientation) .  Thus  the  general  photogram- 
metric  problem  is  essentially  an  interpolation  procedure  by  which  systematic 
errors  may  be  eliminated  and  the  propagation  of  residual  errors  in  the  final 
triangulation  procedure  decreased.  Depending  upon  the  arrangement  of  the  triangu¬ 
lating  pencils  of  rays  (normal  case  or  oblique  case)  and  upon  the  type  of  photo- 
grammetric  instrumentation  used  (aerial  cameras  or  phototheodolites),  the  number 
of  independently  given  orientation  elements  will  vary.  However,  this  fact  does 
not  eliminate  the  basic  necessity  of  computing  the  remaining  orientation  elements 
by  a  least  squares  adjustment  as  stated  above.  Thus,  it  appears  to  be  unnecessary 
to  distinguish  between  different  types  of  photogramnetry  because  the  methods  of 
atereophotogrammetry  applied  to  all  measuring  problems,  e.g.,  in  geodetic  and 
industrial  work,  are  based  on  essentially  the  same  geometrical  and  physical  prin¬ 
ciples. 

A  general  theory  for  the  error  propagation  deals  with  the  mean  error  of  an 
observation  of  unit  weight,  with  the  errors  cf  the  elements  of  orientation  and 
their  propagation  into  the  mean  errors  of  the  spatial  coordinates  of  the  model. 


II.  GENERAL  REMARKS 

The  concept  of  measuring  should  always  include  the  principle  of  overdeter - 
mlnation.  Therefore,  a  final  determination  of  the  unknown  parameters  should  be  the 
result  of  an  adjustment.  The  significance  of  such  an  answer  is  determined  by  the 
computed  estimates  of  the  precision.  In  general,  it  is  possible  to  obtain  approxi¬ 
mate  values  for  the  unknown  parameters  without  difficulty.  Consequently,  in  the 
process  of  adjustment,  corrections  to  the  approximate  values  must  be  computed  to 
render  the  final  result.  The  process  of  a  rigorous  adjustment  requires  that  the 
weighted  sum  of  the  squares  of  the  residuals  of  the  original  measurements  be  a  min 
imam.  If  the  measuring  errors  are  normally  distributed,  the  result  represents  the 
most,  probable  values  for  the  unknown  parameters.  The  precision  of  the  result  is 
computed  from  the  weight  coefficients  obtained  from  the  adjustment  and  the  mean  er 
ror  of  unit  weight,  which  way  be  either  independently  determined  or  computed  from 
the  sum  of  the  squares,  cf  the  residuals .  The  specific  least  squares  adjustment  is 
given  by  rigorous  roathematical  expressions  which  connect  the  original  measurements 
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with  the  unknown  parameters.  If  these  functions  are  not  linear,  it  is  convenient  to 
linearize  them  by  applying  the  Taylor  Beries,  neglecting  second  and  higher  order 
terms .  In  such  cases  it  may  he  necessary  to  compute  a  series  of  corrections  by  on 
iterative  process  in  order  to  achieve  a  desired  accuracy  in  the  final  result. 

In  an  analytical  treatment  of  any  triangulation  problem  -  including  photo - 
grammetric  problems  -  the  above  mentioned  procedure  must  be  considered  as  a  part 
of  a  least  squares  adjustment. 

Before  analyzing  out  problem  in  detail,  let  ua  state  the  objective.  The  sole 
purpose  in  topographic  photogrammetry,  and  in  numerous  non -topographic  applications 
of  photogrammetric  measuring  methods,  is  the  reconstruction  of  the  object  photo¬ 
graphed.  In  accordance  with  such  an  objective  the  efforts  in  research  and  practi¬ 
cal  application  are  often  limited  to  problems  of  production  and  quality  of  the 
stereoscopic  model  and  of  its  evaluation.  In  such  cases,-  a  limited  amount  of  empha¬ 
sis  is  put  on  the  geometry  and  the  physical  conditions  existing  during  the  process 
of  taking  the  photographs.  In  various  non-topographic  applications,  however,  photo¬ 
graphs  are  taken  either  to  calibrate  certain  physical  parameters  of  the  cameras 
themselves  or  to  determine  absolute  values  of  the  orientation  elements-*-.  Typical 
examples  are:  a)  the  determination  of  position  and  attitude  of  an  airborne  camera 
from  photographed  ground  control  points  for  the  purpose  of  trajectory  determina¬ 
tion  of  airplanes  or  grided  misriltas,2  and  b)  the  calibration  and  orientation  of 
photographs  taken  with  ground  baaed  cameras  for  the  purpose  of  providing  a  refer¬ 
ence  datum,  where  Buch  cameras  are  mounted  rigidly  in  relation  to  other  measuring 
equipment.  In  addition,  a  rigorous  mathematical  treatment  of  the  problem  of  tri- 
angulation  is  necessary  for  a  study  of  the  propagation  of  random  and  systematic 
errors  in  both  single  models  and  extensions. 

In  general,  an  approximate  solution,  though  useful  for  a  specific  case,  will 
rot  necessarily  be  adaptable  to  other  cases.  However,  a  general  and  rigorous  ana¬ 
lytical  treatment  of  the  simultaneous  orientation  of  two  photogrammetric  cameraa 
will  be  applicable  to  any  measuring  method  based  on  the  Btereophotogrammetric 
principle.  Approximate  solutions  which  may  he  desirable  for  economic  reasons  in 
special  cases  can  be  derived  from  a  general  solution  by  introducing  certain 
assumptions  peculiar  to  such  cases.  Therefore  an  analytical  treatment  should  be 
sufficiently  general  as  to  the  least  squares  procedure  and  the  code  for  electronic 
computers  to  permit,  the  use  of  the  method  for  a  vide  range  of  stereophotOgrammetric 
measuring  problems. 


III.  THE  GEOMETRICAL  ANALYSIS 

In  making  a  rigorous  geometrical  analysis  of  the  simultaneous  orientation  of 
two  photogrammetric  cameras,  we  shall  first  consider  the  method  of  deriving  the 
orientation  used  in  any  optical -mechanical  stereoscopic  restitution  instrument . 

The  spatial  position  and  orientation  of  a  camera  with  known  optical  characteristics 
are  defined  by  six  parameters,  namely,  three  translations  and  three  rotations. 
Consequently,  for  the  two  cameras  under  consideration  we  are  dealing  with  12  un¬ 
known  parameters . 

Two  bundles  of  projective  rays  can  be  oriented  with  respect  to  each  other 

1  cf .  [ 8] .  ( References  at  the  end  of  the  paper) 

2  cf.  [9]. 


from  measurements  made  on  the  photographs  without  reference  to  absolute  control 
data.  Five  pairs. of  corresponding  rays  are  necessary  and  sufficient  to  establish 
this  relative  orientation  which  is  equivalent  to  the  construction  of  a  true  model 
of  the  object  photographed  by  the  two  cameras.  In  the  process- of  determining  the 
remaining  seven  unknowns  we  must  find  the  scale  and  the  absolute  orientation  of 
the  model.  The  latter  establishes  the  model  with  respect  to  a  given  control 
system  by  three  additional  rotations  and  translations  which  are  equivalent  to  a 
coordinate  transformation. 

The  separation  of  the  orientation  problem  into  the  two v fundamental  processes 
of  relative  and  absolute  orientation  has  been  predominant  in  practical  photo- 
grammetiy  and,  consequently,  has  greatly  influenced  the  related  error  theory. 

The  separation  of  the  two  processes  of  orientation  has  been  strictly  based  upon  the 
practical  methods  of  restitution.  From  the  standpoint  of  both  practice  and  error 
theory,  this  separation  has  its  disadvantages.  The  model  obtained  from  the 
relative  orientation  -  even  if  it  were  flawless  -  is  of  interest  in  '' very 
special  cases  only.  The  relative  model  usually  has  to  be  transformed  into  an 
absolute  one.  During  this  process  absolute  control  in  excess  is  often  used  to 
improve  the  preliminary  relative  orientation  (compensation  for  model  deformations) 
as  well  as  to  minimize  the  influence  of  existing  tensions  within  the  given 
control.  In  any  case  the  absolute  control  must  be  considered  to  be  of  greater 
inport ance  than  the  relative  information. 

In  an  analytical  treatment  the  absolute  orientations  of  the  two  projective 
bundles  of  rays  under  consideration  are  functions  of  the  positions  of  photographed 
images  on  the  plates.  This  conception  automatically  includes  the  relative  as 
well  as  the  absolute  orientation  of  the  model..  In  addition,  the  sum  of  the 
squares  of  the  weighted  residuals  of  the  original  coordinate  measurements  on  the 
photographs  must  be  minimized.  Therefore,  the  orientation  problem  may  be  expressed 
by  formulas  explicit  in  terms  of  the  measured  plate  coordinates  and  preferably 
so  arranged  that  a  minimum  number  of  measurements  and  residuals  appears  in  each 
of  the  observational  equations. 

It  is  of  further  interest  to  note  that  in  an  analytical  treatment  the  method 
of  observation  -  monocular  or  stereoscopic  -  does  not  influence  the  formulas 
..^pressing  the  rigorous  geometry.  In  other  words,  monocular ly  and/or  stereoscopical- 
ly  observed  coordinates  may  be  introduced  into  the  computation  so  long  as  the 
condition  i3  satisfied  that  the  sum  of  the  squares  of  the  weighted  residuals  of  the 
original  measurements  is  minimized.  Consequently.,  an  analytical  method  may  use 
absolute  control  data  outside  the  field  common  to  both  photographs.  Thus  an 
extended  basal  area  is  obtained  for  the  double  point  resection  in  3pace,  a  fact 
which  will  make  it  possible  to  use  the  favorable  geometry  which  is  characteristic 
of  the  single  spatial  resection-^ . 

In  order  to  obtain  a  general  solution  for  our  problem,  the  least  squares 
adjustment  must  be  suitable  for  all  numbers  and  types  of  reference  points  and 
their  corresponding  plate  measurements  -  provided  that  the  information  is 
sufficient  for  a  unique  solution  and  satisfies  the  principle  of  the  relative 
orientation.  The  three  types  of  reference  points  are*.  1)  absolute  control 
points  given  by  their  spatial  coordinates,  2)  partially  absolute  control  points 
given  either  by  their  Z  coordinates  or  cy  their  X  and  Y  coordinates,  and  3) 
relative  control  points.  Points  of  the  first  type  are  not  necessarily  restricted 
to  the  area  common  to  both  photographs  but  the  two  latter  types  of  points  must 
lie  within  this  area. 


The  basic  geometry  necessary  to  meet  the  above  mentioned  requirements  is 
simple.  Assuming  a  unique  solution,  we  see  from  Figure  1  that  the  condition  for, 
correct  orientation  with  respect  to  an  absolute  control  point  R>  given  by  X,  Y,  Z, 
is  equivalent  to  the  condition  that  the  center  of 'projection  0,  the  image  point  r 
and  the  control  point  R  are  colllnear.  In  case  the  absolute  control  point  la  located 
within  the  area  common  to  both  photograph*,  the  condition  of  colli nearity  Is  valid 
for  both  the  left  and  right  photographs.  In  such  a  case,  the  condition  of  inter¬ 
section  of  two  corresponding  rays  -  namely  at  point  R  -  1*  automatically  satisf  ied. 
The  fact  that  the  left  and  right  cameras  are  fully  independent  in  regard  to  the 
condition  of  colllnearlty  simplifies  the  least  squares  adjustment.  Furthermore, 
this  fact  is  interesting  from  the  standpoint  of  error  theory.  A  similar  geometri¬ 
cal  condition  may  he  obtained  for  relative  control  points .  Also  from  Figure  1  we 
see  that  the  condition  that  two  corresponding  rays  intersect  is  equivalent  to  the 
condition  that  the  two  center#  of  projection  O'  and  0",  together  with  two  corre¬ 
sponding  image  points  r '  and  r" ,  are  ooplanar.  In  deriving  the  necessary  formulas 
we  shall  sea  that  the  two  geometrical  conditions  mentioned  above  are  Adequate  for 
the  aolution  of  our  problem. 


From  Figure  1  we  obtain: 

n i  .  . .  |Mi 
St  —  A 

"V  +  ( l) 

and  R"  -  u"r” 

where  n*  and  pi"  sue  scale  factors.  The  prime  and  double  prime  lndioee  refer  to 
left  and  right  stations,  respectively. 


Correspondingly  i 

X  •  XJ  ♦  n'm*  •  XJ  ♦ 

X  •  YJ  ♦  p»v»  •  I«  ♦  |4»v«  (8) 

Z  ■  Z1  ♦  uhs*  •  Z"  ♦  n^w" 
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Snch  of  the  triplets  of  formulae  (2)  is  the  analytical  expreaalon  for  the  condition 
that  the  points  O',  r'.  |R  and  0",  r",  R,  respectively,  lie  on  straight  lines.  Be- 
oauae  the  point  R  la  ooamoo  to  both  triplets,  these  lines  satisfy  in  addition  the 
condition  of  Intersection  at  the  point  R. 


Frcm  fomila  (2)  follows: 


X.Ia.  (!) 


V  -  ^  and  t*  -  where  X  -  Tn  -  (Y)  (3) 

w«  -  Z-Z  •  (Z) 

u  p"  o 


Free*  Figure  2  «ie  road 

F  ■  iu  >  iv  ♦  kw  •  "£x  ♦  3y  ♦.lie 
R  -  i(X)  ♦  J(T)  ♦  k(Z) 


iu) 


9 


;  V 

Further  vs  obtain  r  It  -  c  and  by  formulas  (l)  ^  r  k  *  R  k  and,  therefore, 

R  it 

With  the  transformation  matrix  of  tha  two  veotor  triplet* 


and  fraa  formulas  (U),  v«  obtain 

9  .  ^  ■  ■■ 

m  (X)  tin  a  ooa  «a  ♦  (T)  ain  m  ■»  (Z)  ooa  a  ooa  m 
**  "  o 

““  \\ 
u  ■  x  ooa  a  «y  ain  a  ain  to  ♦  o  alb  «  ooa  o» 

t  ■  y  ooa  «  ♦  o  tin  to 

w  -  -x  ain  o  -  y  ooa  a  ain  *  ♦  o  ooa  a  oos  » 

In  fonula  (7)  «*  hams  from  Figura  J 

5  •  -(x  -  x  >  osai  -  {y  «  y_)  *1»k 

v  P 

7  •  »(x  -  Xp)  Bin  it  ♦  <7  -  7p)  ooa* 

V*  now  lr.troduoa  auxiliaries  which  era  used  throughout  tha  report* 


(5) 


(6) 


(7) 


(8) 


•  -ooa  o  ooa  ♦  sin  a  ain  «  ala 

B.  ■  -ooa  «  ein< 

1  i; 

c  s  Sin  £  *6!  s  —  km  a  art  n  a  tri*  v 

1  \  MM—  —  - -  V 

•  sin  a  cos  to 


*1- 


sin  to 


F  •  cos  a  coo  to 
1  ,  ./ 

0^  «  sin  a  „ 

H-  **  sin  * 
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Aj  »  -oo#  a  sin  <  -  sin  a  ain  m  ooa  ^ 


ooa  <*  ooa  x. 


C.  •  ain  o  ain  x  -  ooa  a  ain  »  ooa  x, 


u2 

*2 


ain  a  ain  •» 
ooa  • 


( 9 ) 


?2  ■  cos  a  sin  to 
C>2  •  cos  a 
H?  »  cos  x. 

-  \ 


By  substituting  (9)  into  (6) 

U)\  ♦  (!)%  4  (Z)?1 


and  by  substituting  (8)  into  (7)>  using  the  expressions  (9) 
a  •  (x  -  Xp)^  ♦  (y  -  yp)A2  * 

t  »  (*  -  *p)R^  ♦  (y  -  Jp)B2  '*  °^i 


w  •  (x  - 


*p>Ci  4  {y  *  7P)C2  *  ®Y 


Substituting  foreulas  (10)  and  (11)  into  (2)  end  (3)#  and  using 


X  -1122  f.  T  snd 

WO 

we  hen  for  the  left  end  right  stations* 


T  miih.  ♦! 


T  (a-z«)  [<x,"*b)Al *  v  0>Dil  Tt 

(x»-*£)CJ  4  (y*-7j)C$  ♦  o»FJ  *  0 

(Z-ZJJ)  [(x«-xJ)A«  4  (y«-yjj)A|  ♦  o»Dj] 

- - —  ♦  X» 

(x"-xJ)CJ  4  (y«~yJ)C|  ♦  o»ri[  0 

( Z-Z£ )  [  (x ' -x^)Bj  4  (y’-y^Bj  4  q»E£] 

(x'-x£)C{  4  (y’-y^)CJJ  40>yj  0 


(Z  -  2")  [(x"-tJ)B»  4 


4  0"«l 


whore 


(x"~xJ)C«  4  (y"-yJ)Cg  4  c"F|  e 

. .  «■  ■>  *  <z-2;>°il 

x  -  -  q?  4  Xp 

c«  [Of-X')A£  4  (T-Y')BJ  4  (Z-ZMCll 

7  o'  yp 

-  (X-X^)r,|  -  (Y-fJ)Ej  -  (Z-Z')F| 
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where 

q«  -  (X-X«)D«  ♦  (T-T*)SJ  ♦  (Z«Z£)F| 

k  more  detailed  derivation  of  the  formulas  (13)  and  (lU)  any  be  found  inLGi)^* 

The  coordinates  of  R  nay  be  only  partially  given,  either  with  its  elevation 
i  or  with  its  location  X  and  Y.  ■  V*  obtain  fro*  the  three  equalities  in 
formula  (2),  in  case  Z  is  given 

(35) 


(16) 


In  formulas  (15)  and  (16)  we  have  Introduced  the  base  line  ooiqpGmnts 


b  -  X"  -  I* 

x  o  o 


b  ■  T*  •*  Y' 
y  o  o 


b  •  Z»  -  2‘ 

*  o  o 


(17) 


Formulas  (15)  are  obtair*d  by  equalising  the  X  and  Y  coordinates  of  point  S 
from  the  left  and  right  stations.  Thus,  for  a  given  value  of  Z,  the  condition 
of  intersection  of  th  two  corresponding  rays  is  satisfied  Similarly,  fonailas 
(16)  are  obtained  by  equalizing  the  Z  coordinate  of  point  R  froa  the  left  and 
right  station.  If  we  introduce  j>he  given  values  for  the  X  and  T  coordinates, 
the  condition  of  intersection  is  again  satisfied.  Bach  equation  in  for«rlas  (15) 
and  (16)  raay  re  written  a<*  a  conditional  equation.  From  formula  (15)  We  obtains 


(Z)',uH¥'  -  (7)  *n'v*  +■  b  v'w"  ■  0 

and  *  (3 

(Z)»v"w'  -  (2) 'vV  ♦  b  w'w"  •  0 

Fron  fomu las  (16)  wa  derive! 

(X^u^w1  -  (x)*u'w*  -  b^tt'a*  •  0 

and  * .  (3 

(l)‘v*W  -  (Y)»v»w»  -  bgv*y"  -  0 

Finally,  we  shall  consider  the  case  in  whioh  the  absolute  position  of 
paint  R  is  unknown  and  only  the  condition  of  intersection  is  given*  Fros  the 
equality  in  formula  (1$)  we  obtain! 

b  (v'v"  -  v***)  ♦  b  (u*w'  -  S’w*)  .♦  bg(u<v*  -  w"v»|  •  0  (2 

"  5 

This  equation  may  be  directly  obtalnod  by  a  different  approaoh*  ..  The  condi¬ 
tion  of  int*raaotion  for  any  two  corresponding  rays  is  satisfied  if  the  two 
rays  lie  in  a  plane  (baeal  plane)  with  ihe  base  line  O'O".  (See  Figure  1). 
This  condition  is  equivalent  to  the  coixiition  that  the  four  points  O',  0*,  r* 
and  r*  are  coplanar*  The  coplanarity  of  these  four  points  is  satisfied  if 
the  determinant  formed  by  the  spatial  coordinates  vanishes.  With  respeot  to 
the  (X),  (Y),  (2)  -  system  we  have!  » 


b^  ♦  un 


by  ♦  V* 


bg  ♦ 


This  determinant  is  identical  with  the  conditional  equation  (20)*  Foaula 
(20)  represent*  the  basic  expression  for  the  relative  orientation  for  the 
general  case.  If  both  b  and  b  are  equal  to  aero*  ye  have i 


or  from  formulae  (12) 


For  the  case  in  which  b  •  0,  the  elevations  of  the  stations  2'  and  Z*  are 

»  oo 

equal,  and  consequently  (2)'  -  (2)*.  Thus  we  may  write  ( T) •  ■  (l)“,  which  is 

the  traditional  presentation  for  tie  condition  of  relative  orientation  as 

given  by  v.  Oruber.6  The  conditicn  of  equality  for  the  (T)  coordinates  is 

used  in  the  process  of  relative  orientation  of  independent  pairs  of  photographs 

in  the  optical-mechanical  restitution  equipment.  Actually,  the  observation 

of  any  (I)  coordinate  makes  it  necesoary  to  introduce  a  corresponding  (Z) 

coordinate.  Consequently,  tuo  {.-iOCbSS  of  relative  orientation  is  sctus-lly  tssed 

on  formula  (22)  wr.:  oh  may  be  written  esi 


-f*  N 

:f .  [3]  pp.  27  -  2y 
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H 


t.'.n  X  •  ■  tan  X*  (23) 

X  is  the  ipatial  -^ngle  formed  by  an  arbitrarily  chosen  reference  plane 
containing  the  base  line  and  a  second  plane  containing  the  base  line  and 
the  corresponding  target  ray.  Thus,  X  is  the  polar  bearing  of  K.  Equation 
(23)  states  that  the  two  polar  bearings  of  R  as  obtained  from  the  left  and 
right  stations  must  be  equal,  a  requirement  which  is  identical  with  the  above 
mentioned  general  condition  that  the  two  corresponding  rays,  with  the  base 
line,  lie  in  a  plane.  In  this  connection,  it  may  be  mentioned  that  e.g.  in 
Finsterwalder'e'  formula  (10),  given  as  the  basic  formula  for  the  relative 
orientation,  the  first  quotient  equals  tan  X"  and  the  second  quotient  equals 
tan  X'.  Thus,  his  formula  represents  the  special  case  of  the  relative 
orientation  as  expressed  by  our  formulas  (zjj. 

Summarising  the  results  obtained  so  f  ar,  we  see  thati 

a)  Each  absolute  control  point,  defined  by  its  X,7,2  coordinates,  gives 
rise  to  two  independent  equations  for  each  nodal  point  position  from  whioh  the 
oontrol  point  is  photographed*  Thus,  if  the  absolute  control  point  lies  in 
the  area  common  to  both  photographs,  we  obtain  four  independent  equations 

and  in  addition  satisfy  the  condition  that  the  two  corresponding  rays  intersect* 
(For nu las  (ll»)). 

b)  Each  partially  absolute  oontrol  point,  given  by  either  its  2 
coordinate  or  its  X  and  Y  coordinates,  gives  rise  to  two  Independent  aquations 
whioh  satisfy  the  condition  of  intersection  of  two  corresponding  rays* 

(Formulas  (10)  or  (19),  respectively). 

o)  Eaoh  relative  control  point  gives  rlss  to  ont  independent  equation 
and  represents  the  condition  that  two  corresponding  ray?  intersect.  (Femaiia 
(20)). 


From  the  results  given  In  (e)  to  (c)  above  we  conclude  that  for  a  unique 
solution  it  is  nocessery  and  sufficient  to  have,  for  exanple,  for  the 
determination  of  twelve  unknowns!  a)  two  absolute  oontrol  points  whioh  are 
situated  In  the  area  common  to  both  photographs  giving  riso  to  eight 
equations  and  satisfying  two  conditions  of  intersection]  b)  cos  partially 
absolute  control  point  giving  rlss  to  two  additional  equations  and  increasing 
the  riunbor  of  enforced  intersections  to  three)  o)  In  addition  two  relative 
oontrol  points  increasing  the  nuntoer  of  equations  to  the  necessary  twelve  and 
establishing  the  fourth  and  fifth  intersections  of  corresponding  pairs  of  rajt* 

These  requirements  are  identical  with  those  necessary  and  sufficient  to 
solve  the  problem  in  the  conventional  way,  widch  consists  of  dividing  the 
orientation  prvblum  into  iha  phases  of  relative  and  absolute  orientation  as 
stated  at  the  beginning  of  this  chapter. 

In  applied  aerial  phot  ograrcre  try  the  difference  in  flying  height  between 
two  successive  photographs  -  the  bane  cor^ponePt  bs  -  is  often  determined 

independently  >  y  Statoscopo  registration,  and  experience  has  shown  that  It  is 
s  ore  titles  desirable  to  enforce  U.;r.  measurement  in  the  process  of  orientation. 
Analytical  ,y  speaking,  t(  ere  exi-its  in  such  a  case  an  additional  conditional 

nf  tl’,‘  fC’JT'l 

'll 


Z"  -  Z '  -  b  -  0 
0  0  8 

Because  of  the  simplicity  of  this  formula  it  seems  advisable  to  write 


1*  -  Z •  ♦  b 
0  o  s 


(21*,) 


and  substitute  this  expression  into  the  formulas^  (1U),  (18),  (19),  and  (20). 


Similarly,  the  same  procedure  may  be  followed  for  the  b^  and  by  components* 
The  formulas  introducing  the  unknowns  bx,  by  and  b#  are  given  below, 

. .  <”  [W'm  ♦  « <f>,ci3.  . 

X  m  *p 

c«  [  (X)'A'  ♦  (Y)'BA  ♦  (Z) *Ci] 

r  -  — - - f  -ji— 2 — 

q'  -  (X)»D»  ♦  (Y)»BJ  ♦  (Z)'F£ 

o» {((!)•  bj  A$  ♦  [<Y)‘  -  b  ]  Bl[  ♦  [(Z)'  -  bj  Cj} 

V 

o»  {[<*)'  -  bx]  A|  ♦  [(Y)«  -  by]  Bg  ♦  [(g) «  -  bj  eg} 


where 


and 


<»") 


*n 


where 

q"  -  [(X)' 

-b*] 

W[  ♦  [(!)• 

(Z)i(u**w‘  - 

W'w") 

-  b_u*w'  ♦ 

and 

m 

(Z)'(v»w‘  - 

Vw") 

m  b#V*W»  ♦ 

(X)'(u*w«  - 

u»w") 

♦  bxu'w"  - 

(Y)  - 

v*w") 

♦  byV’w*  - 

Formula  (?Q)  remains  unchanged. 


- »,]  *s  *  [(2)'  -yi  q 


b  w»w"  •  0  / 

(lfi*) 

b  w'w»  ■  0  S 

7 

b.u'u*  ■  0 

)  ,  . 

s 

(19*) 

b^T'V*  ■  0 

In  order  to  enforce  ary  or  all  base  line  conponents  in  the  orientation 
process  the  corresponding  unknown  parameter  corrections  need  only  to  be 
eliminated  from  the  confutation. 


V 
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IV,  THE  LEAST  SQUARES  ADJUSTMENT 


A.  The  process  of  orientation 

1)  The  least  squares  solution 

Formulas  (li*),  (18),  (19)  and  (20)  express  the  functional  relationship 
between  tile  pie  te  uuOruineted  eiw  the  given  control  data  for  ths  prOCOOS  Of  ' 
orienting  two  photogranmetric  cameras  simultaneously 4  Consequently,  these 
formulas  may  be  considered  as  the  basis  for  a  rigorous  least  squares  adjust¬ 
ment.  In  formulas  (18),  (19)  and  (20)  more  than  one  observation  is  involved, 
leading  to  more  than  one  residual.  Moreover  these  observations  together 
with  the  unknown  parameters  of  the  solution  must  satisfy  certain  conditions. 
Therefore  we  must  deal  with  the  general  problem  of  a  least  squares  adjustment* 
A  complete  treatment  of  this  problem  was  first  prepared  by  He  line  rt8. 

The  form  of  the  linearised  conditional  equations  for  the  different  types 
of  oontrol  points  is  shown  in  the  observational  equations  (2$).  The  actual 
values  of  the”  coefficients  will  be  dealt  with  in  a  later  chapter?. 

Formulas  (29)  show  that  for  absolute  control  points  only  one  observation 
appears  in  eaoh  equation  and  no  observation  appears  in  more  than  one  observa¬ 
tional  equation.  Bence,  the  coefficient  matrix  of  the  residuals  for  these 
equations  is  diagonal  and.  In  fact,  is  equal  to  the  unit  matrix.  Furthermore, 
the  conditional  equations  for  each  station  have  only  those  coefficients  which 
correspond  to  ths  elements  of  orientation  of  their  respective  stations.  As 
already  mentioned  on  page  10,  for  a  case  in  vhioh  no  partially  absolute,  and  no 
relative  oontrol  points  are  present,  the  stereo-orientation  will  reduoe  to  '' 
two  Independent  oamera  orientations* 

In  tha  equations  arising  from  partially  absolute  control  points,  four 
observations  appear  in  each  equation.  However,  the  same  observations  appear 
in  ohly  two  equations.  The  coefficient  matrix  of  the  residuals  Is  therefore 
a  diagonally  arranged  sequence  of  non-overlapping  2x1*  submatrices. 

Finally,  in  the  equations  arising  from  relative  control  points,  four 
observations  appear  in  eaoh  equation  and  none  of  these  appears  in  any  other 
equation*  Thus,  the  coefficient  matrix  of  the  residuals  is  a  diagonally 
arranged  sequence  of  non- overlapping  1x1*  subnmtrices*  ^ 

For  the  partially  absolute  and  relati/e  control  points  the  elements  of 
orientation  for  both  stations  appear  in  the  corresponding  conditional 
equations  and,  consequently,  the  coefficient  matrix  of  the  parameters  is. 
filled. 


8  cf.  ’.5]  pp.  21?  -  222 
^  cf.  pp.  29  -  31* 
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Assuming  the  observations  to  be  independent  and  normally  distributed,  the 
most  probable  values  are  obtained  by  mini m± ting  vTP  V  ,  where  P  denotes  the 
weight  natrix 


P  - 


p1  0 

0  p2 


0 

0 


V. 


(28) 


h-  \ 

oS.  I 


where  the  p's  ere  the  waists  of  the  obeerratl^as. ))  The  introduction  of  weight¬ 
ing  factors  nay:  become  necessary  to  express  vaflous  degrees  of  preoleion  with 
which  the  original  observations  have  been  Obtained  due  to  the  Method  of  manuring 
or  due  to  a  varying  image  quality  caused  s.g.,  by  loss  of  definition  towards 
the  etbrae  of  \$ha  photograph.  In  oase  the  observations  all  have  the  earn 
weight,  it  is  convenient  to  oonsider  P  as  the  unit  aatrix* 

A  direct  solution  of  our  problem  is  obtained,  as  Helasrt  has  shown,  by 
Minimising  tbs  following  funotiom 

¥tPv-8  fcT  (  Av*BA-L  )  (29) 

where h  denotes  a  vector  of  unknown  Lagrange  nultipliara  or  oar relates. 


Setting  each  of  the  differential  quotients  for  Vj,  Vg, .....Vg,  Ap  Ag, 

•  ••••  Au  equal  to  taro,  wa  obtain  (n  ♦  u)  aquations  which,  together  with  the 
previously  mentioned  r  conditional  equations  (formulas  (26)),  are  neoesnaiy  and 
eufflQlant  for  the  determination  of  the  n  reelduala  (v),  the  u  unknown 
parameter  corrections  (A)  end  the  r  oorrelatee  (k).  The  differentiation 
for  v^,  vg, .... »v  gives t 


Pw  -  AV  o 

or  V-  P  '1ATk 

(30) 

The  differentiation  for  ^  p  A  g, . 

....  gives* 

0Tk-  0 

(31) 

Substituting  formulas  (30)  into  (2?)  and  using  formulas  (31),  w#  obtain  the 
(r  *  u)  normal  equations: 


AP*1  ATk*BA- L  -  0 
BTk  -  0 

Because  the  square  Matrix  AP'*AT  la  non-singular  in  our  problem,  we  say 
solve  for  k  in  torma  of  A  . 


(32) 


k-  -  (  AP*lAT  )_1(B  A  -  L  ) 

21 


(33) 


Substituting  (33)  into  (31)  we  obtain  the  final  system  at  notttafck  equations  for 
the  parameter  corrections  2a  * 

[bt(  A  A  -  Bt(AP“1At  r1  L  -  0  Wi 

The  mean  error  of  an  bbservation  of  unit  weight  is 


The  feasibility  of  this  solution  depends  upon  the  fact  that  the  r  x  r  matrix, 
AP'V  ,  is  essentially  diagonal* 

If  the  A  matrix  can  be  partitioned  into  diagonally  arranged  A^  sub-  .. 

matrices,  as  indicated  for  our  problem  with  formulas  (26),  it  oan  be  shown  that 
the  equation  (3U)  may  be  written  as 

L  f*T  ( A  P  '1a'  )_1  ■]  1 A  [>’  (  A  VV  )*1l]  i  -  0  06) 

where  the  index  i  refers  to  the  elements  of  the  1th  group*  Thus  it  is  obvious 
that  for  computational  convenience  ths  final  sat  of  normal  aquations  may  bs 
obtained  by  adding  separata  sets  of  partial  normal  equations  formed  for  eaoh 
of  the  five  grams.  The  sat  of  partial  normal  equations  for  the  1™  group  Is, 
according  to  (36) i 

[aT(  A  P“V  )-1  i]  i  A"  [bt(  A  P-1Ar  )mlL  ]  t  -  0  (37) 

m  analogy,  it  is  possible  to  subpartition  the  matrides  in  (37)  into  units 
corresponding  to  the  Individual  condition  equations  ,’for  a* '  specific  group, 
provided  the  submatrix  A  ^  is  diagonal.  Then  the  set  of  partial  normal 

equations  (37)  may  be  written  as 

^  ■  .  u 

y  • 

[  bt(ap-1aT)-1b]i3  a-  jr  [bt(ap-V)-1l]1j-o  <») 

where  the  index  J  refers  to  the  elements  of  the  line  in  the  group* 

Formula  (36)  shows  that  a  not  of  partial  aornwl  equations  for  a  specific  group 
(formulas  (37) )  oan  be  subpartitioned  into  as  many  sub-sets  of  partial  normal 
equations  ae  the  cormsponding  A  ^  matrix  can  be  subpartitioned  into  diagonally 

arranged  A  ^  submatrices.  Therefore  the  formation  of  the  final  normal 

aquation  aysteu  can  be  accomplished  by  adding  step  by  step  the  partial  normal 
equations  confuted  for  each  of  the  subgroups*  There  is  an  advantage  in  using 
this  procedure  on  electronic  conputers,  because  the  amount  of  memory  space 
needed  is  almost  independent  of  the  number  of  points  carried  in  the  solution. 

-  ti 

In  order  to  conclude  our  investigation  a  study  shall  now  be  made  of  the 
character! sties  of  the  partial  normal  equations  for  each  of  the  smallest 
possible  subdivisions  within  the  five  groups  of  observational  equations  as  „ 
given  by  formulas  (26),  As  mentioned  above  ihe  matricesA^  and  Ag  are 

_  *  22  ..  ■ .  ....  .  .  . 


unit  matrices.  Therefore,  for  both  groups  for  osoh  individual  conditional 
equation  ve  have 

( *  3  *1  AV'1  ■  *3  <»>  ' 

Similarly  in  group  5  the  A  ^  matrix  Is  a  diagonally  arranged  eequenoe  of 

non-overlapping  1x4  suhaa  trices.  Therefore,  ve  have  for  each  Individual 
conditional  aquation  In  thie  group 


Finally,  va  mat  ooosider  groups  3  and  4  which  are  formed  analogously.  Instead 
of  associating  a  sat  of  partial  noraal  equations  with  each  individual  observa¬ 
tional  aquation  ms  is  possible  in  groups  1,  2  and  5*  such  a  sat  nust  ba 
associated  with  a  pair  of  observational  aquations  baoauaa  tha  A .  natrloaa  for 
both  groups  3  and  U  arm  diagonally  arranged  saquanoas  of  non- overlapping  2x4 
aubmatrlooa.  For  tha  J*b  pair  of  obaarvatlonal  aquations,  s.g.  In  group  3, 
va  hare i 


<s)1(a)2(a:i(a)u 
3j  (b^b^btyb)^ 

L  J  33 


Omitting  tha  derivation  and  using  tha  notation  of  formulae  (25),  the  following 
aat  of  partial  noraal  aquations  can  bo  formed  for  eaoh  pair  of  corresponding 
obaarvatlonal  aquations  in  groups  3  and  4,  raapaotlvaly. 

\\ 
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The  weighting  factors  (g)^  (g)2  and  (g)1>2  *n 


and  analogous!/ 


*)T(*)X> 
P  J 


J±tli 

*n  . 


whore  K. 


i(b)(b] 

\  P 


/(aUl)  C\ 

v  iH-/  j  - 


From  (U2)  it  le  obTioua  that  ^his  set  of  partial  normal  equations 
may  be  further  partitioned  by  forming  seta  of  partial  normal  equations 
which  incorporate  only  the  firei*  or  second  terms  of  eaoh  of  the  coefficients* 

For  an  analytical  treatment  the  eeono^r  of  the  least  squares  solution 
is  an  iaqjortant  factor*  Therefore  the  information  so  far  obtained  Is 
summarised  and  the  least  squares  solution  for  the  process  of  orientation 
is  arranged  in  the  following  way.  By  introducing  pseud o-rasiduals  X  we 
"may  denote  in  ('- 


group 
1  and  2 


3  and  U 


$ 


The  observational  equations  (26)  nay  now  be  written  as  follows 


group 


weighting  factors 


®3  A  ' 

Bu 


(9)  and  (9) 
(9)  (9) 


In  terns  of  fomulas  (2$)  the  observational  equations  are  shown  on  the  next 
page. 

The  pseudo-residuals  ((X),  and  {X}  )  are  linear  combinations  of  the 
original  residuals  v  (formulas  (1»6)  )»  They  have  been 

introduced  solely  for  the  purpose  of  sliqpllfying  the  oonputatlonal  procedure. 
The  system  (1*8)  resembles  e  set  of  observational  equations  for  independent 
indirect  measurements.  In  the  last  Qolumn  the  corresponding  weighting  factors 
ars  recorded.  As  the  first  step  in  the  prooess  of  forming  the  normal  equations 
a  sst  of  partial  normal  aquations  is  computed  for  all  r  observational  equation* 
according  to 


(  eVaJi  -  £  <bVl  \ 


P*  denotes  the  diagonally  arranged  weighting  factors,  where  (g)p  (g)2»  (g)p 
(g)y  and  [g)  are  formed  by  formulas  (1*0)  and  (1*5) •  The  index  3  denotes  the 
elements  |#f  the  i**1  line. 


Now  it  Is  necessary  only  to  take  into  account  the  faot  that  in  groups  3 
and  1*  the  pseudo^tesiduals  (X)  by  pairs  are  linear  combinations  of  the  original 
residuals.  Denoting  the  number  of  pairs  of  such  equations  included  in  the  r 
observational  equations  by  s,  another  set  of  partial  normal  equations  may  be 
formed  according  to 

-  Ij  S\j  **1  j  Bjt)  k  A  *  (  B  ji  PV  Lijlk'°  <») 

If-  | 

where  #  *  (g ^i  J  d*no'tea  the  cross  weighting  faotor  for  a  pair  of  , 

corresponding  observational  equations.  The  submatrices  and  B  . .  are  the 

^  J 

matrices  of  the  coefficients  of  the  unknown  parameters  in  the  pair  of  cbserva- 
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tional  equations  under  consideration.  The  two  rows  of  these  matrices  are 
arranged  according  to  the  ij  or  ji  sequence^ of  the  lines  1  and  j,  respectively. 

The  conformity  of  the  two  sets  of  partial  normal  equations  (1*?)  and  (50) 
contributes  markedly  to  the  economy  of  the  numerical  solution.  The  final 
system  of  normal  equations  is  the  algebraic  sum  of  the  two  sets  of  partial 
normal  equations  and  may  be  symbolised  by 


9y  introducing 
we  obtain 


BT  P8A-BtPL-0 

btpb-n 


8tPL 


fy"1  is  the  matrix  of  the  weight  coefficients  of  the  parameters.  The  mean 
error  of  an  observation  of  unit  weight  is  conymted  from  (35) »  where  the 
numerator  is  obtained  during  the  reduction  of  the  normal  equations.  The 
modernised  Oaussian  algorithm*1  is  suggested  for  this  procedure.  The  individual 
v's  for  the  absolute  control  points  are  then  computed  directly  from  the  first 
two  groups  of  the  observational  equations  (1*8)  by  formulas  (1*6)  as 

« ■  x  ($u) 

From  the  third,  fourth  and  fifth  groups  of  equations  (1*6),  the  (X)  and 
fx}  values  are  obtained.  The  individual  v's  are  then  computed  by  formulas  (30), 
(33)  and  (1*7).  Thus  we  obtain  for  a  partial  oontrol  point,  in  the  case  where 


where  ■  *  (^)^(b)^  *  ^^2®1j2 
%  *  "  (x)2^2  *  ^lgl,2 


Px 

pf1"’"" 

y 

(a)^ki  ♦  (o)^kg 
- - - 

<•>!> 


and  in  the  oase  where  X  and  I  are  given* 

vi  -  ^lkl-  *  *Vlk2 

- p - 

f\  yx 

«  ,  m  ^*X?2kl  4  ^aY^2k2 

y  p* 

•>  .  uk.i 


(ax)3ki  ♦  («z)3k2  where  ki  *  -  {X)X(g)X  +  (X)T(g)IT 


kg  -  -  (X)T(g)T  ♦  (X^Cg)^ 


cf.  [ll] 


n 


Finally,  for  relative  control  points  we  obtain! 


#k 


Ty 

i.k 


where  k  •  -  {l} 


During  the  computations  the  usual  checks  are  made: 


W 


jBpv]-  0 


andv^Pv*  [Ll  •  u] 


where  |ll  *  is  obtained  during  the  process  of  reduoing  the  normal  equations* 

The  final  aha  ok  is  obtained  byusing  the  final  orientation  elements  and 
formulas  (14),  (18).  or  (19)  and  (20),  where  either  the  adjusted  observations 
(1  *  v)  and  u'  ♦  v')  mist  oheok  (formulas  (XU) )  or  where,  by  introducing  the 
corresponding  adjusted  observations  together  with  the  final  values  for  the 
unknowns,  the  corresponding  conditional  equations  must  be  satisfied  (formulas 

(18),  (19),  (20)). 

The  mean  errors  of  the  unknown  parameters  are  computed  by  multiplying  the 
mean  error  of  unit  weight,  whioh  oan  be  obtained  from  formula  (35),  by  the 
square  root  of  tw  corresponding  weight  coefficient  whioh  is  obtained  from 
the  matrix  of  the  weight  coefficients  fi|-l.  1«  If 

If  formulas  (14*),  (18*),  (19*)  and  (20)  are  used  as  the  basis  for  a 
least  Squares  treatment,  the  arrangement  of  the  aero-elements  in  the  matrix 
of  the  coefficients  of  the  dbservational  equations  is  changed*  The  correspond¬ 
ing  system  of  observational  equations  is  shown  on  the  following  page*  (formulas  (4th) 

The  increase  in  zero  elements  in  the  third,  fourth  and  fifth  groups  of  the 
observational  equations  in  formulas  (46*)  in  comparison  to  (48)  is  desirable 
slnoe  the  presence  of  weighting  factors  oausea  the  computation  of  the 
ooeffioienta  of  the  normal  aquations  in  these  three  groups  to  be  more  ooiqplicated 
than  that  for  the  first  two  groups.  Moreover,  the  relative  control  points  are 
usually  more  numerous  than  the  absolute  control  points.  Therefore,  the  system 
(48*)  provides  a  more  economical  solution.  The  main  advantage  of  formulas  (48*), 
however,  is  that  any  independently  obtained  base  line  component  may  be  introduced 
in  the  computations  simply  by  eliminating  the  corresponding  db  from  the  least 
squares  computation* 

2)  The  derivation  of  the  coefficients  of  the  observational  equations 

The  setting  up  of  the  observational  equations  (48)  or  (48*)  requires 
computing  the  coefficients  of  the  matrices  of  the  unknown  parameter  corrections 
and  of  the  residuals.  In  addition,  the  absolute  terras  and  the  weights  must 
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Th«i  coefficients  of  the  obserratiouel  equations  for  the  systaa  of  formulas  (U8)  ere  now: 
for  foraoles  (1U)  for  fomila*  (18) 

A  -  c <[l)-{I]  ♦  0)  <*>!  •  ♦  ©»’)  (J>1  - 


■  .  '  "  -  -  -  ..  «, 
.  .  "  '  jj  - 

In  the  base  line  components  are  introduced  as  unknowns,  the  coefficients 
i'or  the  system  of  formulas  (u6“)  are  the  sams  for  formulas  (ike)  aa  for  formulas 
(lk)  except  that  the  arrangement  has  been  somewhat  changed*  For  fornulas 
(l8a),  (I9a)  and  (20),  however,  a. few  coefficient  changes  occurs 


for  formulas  (18* ) 

for  fornulas  (194)  for  fornulas  (20) 

Cd)1  «  • 

(D)x  - 

(dV-  • 

(E) 2  - 

(F) 1  -  -II 

(66) 

U)T  -  ♦  I 
(F)x  - 

(67)  ^  "  * 

{?}  -  * 

(f)2  -  ♦  I 

(f)t  -  • 

the  absolute  terms  are  computed  with 

■■■  . 

-***  “  ~  V 

-&ly  -  7°  -  ly 

-(a)x  -  II*(Z)»  ♦ 

&j  »• 

-(a)2  -  -(I.(Z)i 

Lx]  «•) 

-(a)x  -  II.(X)»  ♦ 

[4]  «• 

-(&)y  •  -(I*(I)» 

tl]  *•) 

-  -I*b  ♦  II*b  ♦  m*b 

x  y  * 

the  weighting  faoters  are  computed  with  fornulas  (1*0)  and  (1*5)  • 

Formula  (20)  expresses  the  general  case  of  the  relative  orientation*  It 
may  be  of  interest  to  consider  the  coefficients  when  conditions  valid  for  the 
normal  case  are  introduced.  With  a  ■  180°,  to  ■  0°  and  0°,  we  have  tram 

formulas  (9)  end  (11) t 

A1  "  4 

1 

A2-  0 

and  u  •  (x  -  x p) 

'  ^  - 

0 

Bg  •  +1 

v  -  (y  -  yp) 

Cl“ 

0 

C«  ■  o 

w  ■  -  c 

*!- 

0 

D2  -  0 

v<‘.  ...  ■ 

,(V 

*1" 

0 

B2  -  +1 

•• 

F1  *  " 

1 

f2  -  o 

•• 

°1  ‘ 

0 

02  -  -1 

a 

«1  “ 

0 

^  -  +1 

■■ 

(68) 


(69) 


o 


and  from  formulas  (60)  and  (8)t 

>  «„•  0  -7  3  *  -<!t  *  V 

*«-•  r  V*5  r  •  *(y  -  7p) 

-v.^"  V  -0 

The  auxiliaries  given  with  formulas  (6l)  are,  with  b  -  ba  *  0,  x  •  y  -  0, 

«»  •  o*  and  b^  •  unityi  ,  ¥  . 

I  •  c(y*  -  y')  [l]  •  0  M  ■  0 

II  «-c(x"  -x*)  \  [2j  •  -c  DQ  •  +o 

m  -  (x»y"  -  x"y‘)  DG  •  -y*  GO  •  y‘ 

* 

Thus,  the  coefficients  (given  in  formulas  (65))  necessary  to  forn  the  correspond¬ 
ing  observational  equation  arei 


w 

m  y*x* 

U) 

»  -y'x* 

(Bi 

-  -(c2  ♦  y'y") 

w 

■  (c2  ♦  y'y")  •  -  {B} 

(c} 

-  ex' 

m 

-  -ox"  ji  s 

(H) 

— W’i-o 

(qV 

(I) 

-  -(a\2  -  *c 

o 

•  -04  ^  *  ”c 

anl  the  corresponding  observational  equation  is,  using  the  last  line  of  fonsulas 
(25)  and  (U8*),  and  formulas  (69) i 

c( v*  -  v')  ■  -y*x'  6a1  ♦  (c2  ♦  y'y*)  Aco'  -  ex'  Ax.' 

yy  •' 

.  ^y'x"  6  b'  •  (o  *  y'y")  &»*  ♦  ox"  A*."  -  o(y"  -  y') 


Introducing  y"  -  y'  •  p  and  v"  -  v»  •  v  .  we  have  sines  J***  y" 

/  s  /  r 

2 

v  .♦£*!  ao'  ♦  x"  -  x'  &  i 1  *  c(l  ♦  )(A«a"  •')  -  p  (71) 

pc  c  ^  c»s  y 

Thfis  equation  is  identical  e.g.  with  R.  Finsterwalder 's  formula  (13)^  for  the 
observational  equation  in  approximately  normal  photography* 


13 

cf.[U)  pp.  15U,  155. 
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D.  Tho  process  of  triangulation 

¥ 

1)  The  triangulation  as  part  of  the  process  of  orientation. 

Like  any  other  trlangulation  procedure,  the  stereophotogramraetrLc  measuring 
r..'thod  is  mainly  a  means  of  determining  the  coordinates  of  oertain  target  points 
by  intersecting  corresponding  rays.  As  distinguished  from  intersection  photo- 
grammetry,  stereophotogrammctry  does  not  necessarily  separate  the  orientation 
procedure  from  the  triangulation  phase.  In  fact  during  the  simultaneous 
orientation  of  two  camera  positions  a  sufficient  number  of  the  target  points 
to  be  determined  are  included  as  relative  control  paints  in  the  process  of 
orientation  in  order  to  improve  the  crver-all  acouracy.  The  inclusion  is  based  on 
the  fact  that  arty  two  rays  which  created  the  images  on  the  respective  plates 
originated  from  a  common  point  and,  consequently,  must  intersect  again  after  the 
orientation  is  restored.  Thus,  by  increasing  the  angular  field  of  the  lens  the 
internal  accuracy  of  the  orientation  is  increased  and,  equally  iiqportant,  the 
presence  and  physical  nature  of  systematic  errors  may  be  determined.  The 
residuals  of  the  plate  roea&uremenla  are  obtained  luring  the  least  squares 
adjustment  for  all  relative  control  points  which  were  included  in  the  orienta¬ 
tion  procedure  and,  thus,  the  adjusted  observations  x  ■  1^  +  v^  and  y  ■  ly  ♦  v^ 

are  available  for  both  stations.  During  the  final  check  of  the  least  squares 
adjustment  for  the  orientation  procedure,  the  corresponding  u,  v,  w  values  for 
both  stations  must  be  computed  in  order  to  check  the  conditional  equation  (20). 
These  values  oan  be  used  in  order  to  obtain  the  corresponding  X,T,Z  coordinates 
directly  from  formulas  (l£)  and  (16),  By  means  of  the  auxiliaries  from  (61) 
we  have  t 
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From  formulas  (72)  and  formulas  (2)  we  obtain 


and 
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In  order  to  obtain  the  more  exact  numerical  answer,  the  confutation  of  the 
scaling  factors  n*  and  )iM  should  be  done  with  auxiliaries  which  do  not  vanish. 
Assuming  ooventional  geometry,  the  use  of  the  auxiliaries  [2],  [£),  and  II  la 
therefore  suggested.  From  formulas  (2)  We  now  haves 
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The  computation  of  the  mean  errors  of  the  final  coordinates  will  be 
discussed  in  a  separate  chapter* 


2)  The  triangulation  as  an  independent  computational  procedure* 

Occasionally,  during  the  process  of  orientation,  it  Is  not  deeirable  to 
cany  as  relative  control  points  all  the  points  whose  coordinates  are  to  be 
determined.  For  those  points  not  included  m  the  process  of  orientation,  an 
independent  coordinate  determination  is  necessary.  The  positions  of  the  corre¬ 
sponding  rays  ere  determined  by  the  elements  of  orientation  as  obtained  from  a 
least  squares  adjustment  and  the  four  measured  plate  coordinates.  It  is 
obvious  that,  because  of  unavoidable  measuring  errors,  these  rays  will  not 
interseot.  They  muot.be  made  to  intersect  so  that  the  sum  of  the  squares  of 
the  corrections  to  be  applied  to  the  original  plate  measurements  Isa  minimum* 
Such  a  least  squares  adjustment  of  conditioned  observations  may  be  based  on 
formulas  (20).  Since  there  is  only  one  conditional  equation  present,  we  have 
only  one  normal  equation  which  is  used  to  determine  the  correlate  k.  Then  the 
residuals  are  easily  computed  according  to  formulas  ($7).  The  final  oheok  is 
obtained  by  introducing  the  adjusted  observations  into  the  origins!  conditional 
equation  (20).  From  here  the  prooedure  of  confuting  the  final  X,?,Z  coordinates 
by  formulas  (7U)  follows  the  steps  outlined  in  the  preceding  paragraph. 

Tha  conditional  equation  for  the  triangulation  adjustment  is,  using 
formulas  (25),  (61),  (65)  and  (69)» 


tVj  ♦  «2v;  ♦  ♦  W|/y  -  CA)  -  0  (75) 

and  the  corresponding  normal  equation  is 

(76) 


(77) 

.ter,  n 

(78) 

As  in  formulas  (57)  the  resdluals  arei 
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(79) 

We  have  the  check) 


vTPv  ■  +  *  k  (6o) 


The  computation  of  the  mean  errors  of  the  final  coordinates  will  be 
discussed  in  the  next  chapter. 

Vi 


C*  Th<s  determination  cf  the  mean  -erhori  of  the  observations,  of  the  elements 
of  orientation  and  of  the  triangulation  results* 

The  mean  error  of  an  observation  of  unit  weight  denoted  by  m  is  coagwted 
from  formulas  (35)  „ 


V  n  -  u 


VTP V may  be  obtained  directly  from  the  reduction  of  the  normal  equations 
according  to. formula  ( 59)  or  by 'adding  the  squares  of  the  individual  weighted 
v  values,  Thus  the  mean  error  of  an  observation  1  before  adjustmaat  is 


*1 


Sometimes  the  confutations  of  m  directly  from  the  original 


• 

measurements  e.g.,  using  the  differences  of  multiple  observations*  may  lead 
to  a  value  of  greatee  physical  significance*  The  discrepancies  between  the 
m  values  confuted  by  the  different  methods  provide  means  to  investigate  the 
presenoe  of  systematic  errors.  1 

„  *  .  -V  *< 

The  inverse  of  the  matrix  of  the  coefficients  in  the  final  normal  equation 
system  (52)  is  the  matrix  of  the  weighting  factors*  The  diagonal  elements  „ 
are  the  squares  of  the  weighting  factors  of  the  corresponding  unknown  para¬ 
meters*  Multiplying  these  weighting  factors  by  m  gives  the  corresponding 
mean  errors* 

Often  it  will  bf  of  interest  to  know  the  mean  error  of  a  function  of  the 
unknown  parameters,-1-**  We  may  consider  the  function 

r(o«,  6>'...y^,  «V  «%  ...  yj)  -  o  (81) 

In  oase  the  function  is  not  linear,  we  apply  the  Taylor  series.  This  gives 
F  »  f  ♦  f,dc» . *  The  f.  values  are  thus  the  partial  differential 

O  l  • 

quotients  of  F  with  respect  to  the  unknown  parameters.  If  we  Introduoe  the 
vector 

(fjf2 . fl8)  -fT  (02) 

the  weight  of  the  function  F  denoted  by  Pp,  together  with  formula  (52),  is 


PF  -  »T  Hmll 

The  mean  error  of  F  denoted  by  is  thus 


v4-  ■*<  > 
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(03) 


(6U) 


By  applying  the  above  procedure,  we  may  determine  the  mean  error  of  any 
function  of  the  orientation  elements.  For  example,  for  t>  the  function  F 
is 


F  -  X"  -  X* 

o  c 


ffl«n 


ih 
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and,  according  to  formula  (81),  we  obtain  with  the  sequence  used  in  formulas  (1*6 ) 


»  +1  All  other  f^  -  values  are  sero# 


The  weight  P.  and  the  mean  error  *.  are  new  computed  by  forailas  (83)  end 
°X  °X 

( 81*),  respectively*  If  the  base  line  eosg>onents  b^,  and  b#  have  been 
considered  as  unknown  parameters  (of*  formulas  (1*8)*  ),  the  terms  yr —  » 

\ 

11  * 

«—  and  'pr-  are  obtained  directly  free*  the  corresponding  diagonal  tanas  of 

*b  o 

y  * 

tbs  matrix  of  ths  weighting  factors* 

Tbs  wan  errors  of  the  triangulated  I,  I,  Z  coordinates  are  denoted  by 
■g,  Bg  and  n^*  The  procedure  of  computing  these  mean  errors  depends  upon 

whether  the  point. was  included  in  the  least  squares  adjustment  of  the  orienta¬ 
tion  proease  as  a  relative  control  point  (see  B/l)  or  whether  the  point  was 
afterwards  determined  by  a  least  squares  adjustment  based  on  four  independent  .... 
plate  measurements  together  with, the  Independently  confuted  elements  of 
orientation  (see  (b/2).  The  latter  case  may  be  treated  by  analogy  to  the  already 
discussed  problem  of  determining  the  mean  error  of  a  function  of  the  unknown 
parameters*  After  the  adjustment  by  formulas  (75-79)  each  of  the  spatial 
coordinates  X,  I,  Z  may  be  expressed  by  formula  (71*)  as  a  function  P  of  the 
elements  of  orientation  end  four  independently  adjusted  plate  measurements* 
Denoting  the  partial  differential  quotients  with  respect  to  the  elements  of 


orientation  and  the  plate  measurements  by  f_,  f«, . ...f, o 
V  respectively,  we  mey  introduce  the  vectors  10 

f-  <'if2 . V 


and 


Fj  and 


-  \oi 

and  pV  (FjFjF^F. ) 

The  propagation  of  the  errors  of  the  elements  of  orientation  into  the 
function  F  is,  with  (81*) 

,  -i  V  i 

•(f)  -  »  (  fnml1  V  -  m  .  Q*(f)  (8< 

The  propagation  of  tha  errors  of  th*  Independently  adjusted  piste  measurements 

is 

V)  •  *  •  “(F)  „  (8S 

where  Q^x  may  be  obtained  in  the  usual  way' during  the  reduction^  the 


(F ) 

anaing 


corresponding  normal  equation*  Since  the  orientation  parameters  are  notN' 
correlated  to  the  plate  measurements  used  in  the  triangulation  phase,  the 
combined  mean  er')ror  of  F  is  confuted  by  the  Gaussian  law  of  Drooaffati  on 


ht)  *  Q(f)> 
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In  approximately  normal  photograph  the  influence  of  the  elements  of 
interior  orientation  on  the  lateral nation  of  X,  7,  2  can  be  negleoted 
because  in  such  cases  the  elements  of  exterior  orientation  compensate  for 
the  effect  of  the  errors  of  the  elements  of  interior  orientation.15 
However,  the  accuracy  of  the  elements  of  interior  orientation  is  important  <t. 
in  interpreting  the  mean  errors  of  the  elements  of  exterior  orientation." 

If  the  point  to  be  triangulated  is  included  as  a  relative  control  point 
in  the  process  of  orientation,  we  must  determine  the  mean  error  of  a  function 
of  the  unknown  parameters  as  well  as  of  the  adjusted  observations.  This 
problem  was  tireated  first  by  Helmart.*' 

For  a  detailed  presentation  of  this  problem  in  natrlx  notation  see  [«Q  • 

>It  follows  that  ve  may  introduce  the  vector 

FT-  (f^»  fo» •••••£■%£»  (FI)' u 

where  the  f  end  F  values  have  the  same  meaning  as  In  forsu3as  (87). 

The  mean  error  of  the  function  F  la  now  given  by 

my  -  V  (92) 

where  R  denotes  the  covariance  matrix  of  the  vector  of  the  adjusted  para- 
asters  and  the  adjusted  observations.  The  numerical  treatment  of  this 
problem  is  quite  cumber  soma. 

The  analytical  method  permits  the  use  of  a  large  number  of  relative  and 
absolute  control  points.  The  matrix  of  the  weighting  coefficients  u“*  (e,g. 
formula  (88))  decreases  as  the  overdetermi nation  increases.  Thus  with  a 
Strong  overdetermination  the  result  obtained  from  formula  (90)  provides  a 
satisfactory  approximation  for  the  determination  of  the  mean  orors  of  ths 
triangulated  coordinates,  making  the  utilisation  of  formula  (92)  unnecessary, 

V.  THE  APPLICATION  OF  THE  ANALYTICAL  METHOD 
TO  THE  PROBLEM  OF  CONTROL  EXTENSION 

The  process  of  strip  triangulation  on  the  restitution  instruments  is  based 
upon  the  technique  of  orienting '’individual  photographs  so  that  the  orientation 
of  a  preceding  photograph  and  the  spatial  position  of  at  least  one  point  of 
the  preceding  model  are  enforced.  This  procedure,  as  well  as  the  separation  of 
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relative  and  absolute  orientation  in  the  restitution  of  an  individual  nodal, 
is  dictated  by  the  physical  properties  of  the  restitution  equipment.  The 
enforoemsnt  of  the  above  mentioned  parameters  in  an  aerial  triangulation  is 
comparable  to  a  procedure  which  breaks  up  the  least  squares  adjustment  of  a 
triangulation  net  into  a  series  of  adjustments  of  individual  configurations, 
enforcing  the  result  of  one  configuration,  with  respect  to  orientation  and 
scale,  into  the  following  one*  However,  the  application  of  such  a  procedure 
does  not  meet  the  error  theoretical  requirements.  In  order  to  carry  the 
orientation  and  scale  over  a  strip  of  overlapping  photographs  it  is  necessary 
to  satisfy  the  condition  that  rays  originating  from  three  consecutive  camera 
Stations  intersect  for  at  least  one  point  located  in  the  area  common  to  three 
photographs.  Analytically,  such  a  poi$t  gives  rise  to  three  independent  condi¬ 
tional  equations.  If' we  denote  three  consecutive  photographs  by  (n-1),  n,  (n+1), 
an  obvious  method  of  setting  up  the  three  conditional  equations  is  the  triple 
application  of  formula  (20)  which  expresses  the  condition  of  intersection  for 
pairs  of  rays  originating  from  the  stations  (n-1)  and  n,  n  and  (n+1),  (ivl)  and 
(n-rl).  However,  this  method  does  not  lead  td  a  result  if  the  three  rays  are 
located  in  or  close  to  a  plane.  This  situation  is  common  in  strip  txlangulation. 
From  Figure  1*  we  see  that  threo  independent  conditional  equations  can  be  ^ 
obtained  by  applying  twice  the  above  mentioned  condition  of  intersection 
expressed  by  formula  (20),  e.g.  for  (n-1)  and  n,  and  n  and  (n+1). 

The  condition  that  tnese  two  intersections  occur  at  a  common  point  is 
equal  to  the  condition  that  the  vector  ia  common  to  the  two  loops 

^  -  b n  •  o/and  Rn  -  -  b^  ■  0,  Thus,  with  formulas  (1)  we 

have  the  condition  that  in  the  first  loop  must  equal  the  corresponding 

expression  in  the  second  loop.  Because  the  vector  ?n  is  identical  In  both 

loops,  an  Independent  conditional  equation  is  obtained  by  equalising  the 
soale  factors  obtained  from  the  two  loops.  Denoting  the  components  of 

the  first  and  second  loops  by  the  indices  1  and  2,  respectively,  we  havat  „ 

(93) 


Formulas  (73)  give  for  each  scale  factor  p  three  quotients  correspond¬ 
ing  to  the  projection  of  the  vector U  into  the  coordinate  planes.  If  we 
arrange  the  ground  reference  in  such  a  way  that  the  X-axis  point* approximate¬ 
ly  in  the  direction  of  the  strip,  it  is  obvious  that  the  terms  Involved  "in  the 
auxiliaries [2j  ,  [5]  and  II,  respectively,  (formulas  (61))  will  always  be 
different  from  sero.  Therefore  we  choose  as  the  third  conditional  equation  usln 
(93),  (73)  and  (61) 

[2]  +[5[)  In2  -  0  •  •:  ;/•  v  (9U) 

The  coefficients  of  the  corresponding  observational  equation  are  obtained 
from  formula  (9h)  as  the  partial  differential  quotients  with  respect  to  the 
unknown  parameters. 
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The  corresponding  value  la  computed  from  (92) 

-{4  •  Mgllj  ♦  [50^2  <**> 

It  should  bo  noted  that  all  the  oonponente  uaad  to  form  the  auxiliaries  given 
by  formula  a  (9$)  are  already  available  from  the  setting  up  of  the  other  tiro 
observational  equations  based  on  formulas  (20). 

An  analytical  treatment  of  the  two  methods  of  extension  shows  the 
difference  between  a  simulation  by  instruments  and  the  least  squares  adjust¬ 
ment.  Assuming  a  strip  of  Approximately  vertical  photographs  flown  with 
2/3  overlap,  the  unique  solution  may  be  considered  first.  The  oontrol  is 
schematically  arranged  as  shown  in  Figure  $  for  n  ■  6  photographs. 
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From  Figure  £  we  obtain  the  foil owing  independent  observational  equations* 


point  no. 

• 

type  of  control 
point  ^ 

number  of 
observational 
equations 

; 

number  of 
conditions  of 
intersection 
satisfied 

formlas 
simlation 
bF  . 

instrument 

describing 

least  scares 
adjustment 
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25 

o 

this  tabu la  t  Ion  shows  that  in  both  cases  we  have  36  observational  equations  for 
the.  determination  of  /.the  noo6  orientation  parameters,  n  being  6  in  our  ease.  „ 
Furthermore, v;five  pairs  of  hay*  within  each  pair  of  consecutive  bundles  must 
Intersect  in  order  to  satisfy  the  rigorous  goomatry.  Thus,  if  n  denotes  the 
number  of  photographs  in  the  strip,  wo? need  (n-1)  •  5  conditions  of  intersection. 
For. a  »*  6,  twenty-five  conditions  of  intersection  must  be  satisfied,as  shown 
in  the  tabulation.  ",  -..••••'T'  .V 


In  Figure  6  the  36  observational  equations  are,  'shown. 


From  Figure  5  we  see  that  five  points  are  common  to  two  photographs  in 
each  of  the  end  models  and  ten  points  are  common  to  three  photographs  in  the 
central  portion  of  the  strip.  •  Consequently,  allowing  two  residuals  for  each 
iftage.  point,  we  havofr(5«2  ♦10*3)  *  80  residuals,  .  In  Figure  6  these  residuals 
&m  repi'eaehted  in  the  cbm; spohdi ng  A  matrix  by  opbn  circles  or  dots,  respect-;, 
ivaly,  depending  upon;  whether  the  corresponding  observational  equation  belongs 
to  the? least  squares  Solution  or  to  a  set  of  equations  simulating  the  extension 
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procedure  as  performed  on  restitution  equipment*  Correspondingly,  the  sequence 

of  the  coefficients  of  the  unknown  parameter  corrections  ia  given  by  thin 

and'hoewy  lines  in  the  B  matrix*  It  ia  obvious  that  the  conventional  sethod 

of  strip  triangulation  is  divided  into  a  sequence  of  independent  extensionai 

thus  allowing  for  the  confutation  of  the  six  unknown  elements  of  orientation 

for  each  consecutive  photograph  separately  (heavy  lines)*  Correspondingly* 

the  A  matrix  is  a  sequence  of  diagonally  arranged  non-overlapping  raw  vectors 

(dots).  However,  such  treatment  will  be  correct  only  in  the  case  of  a  unique 

solution  where  all  residuals  are  equal  to  rero  and,  therefore,  the  numerical 

solution  will  not  be  influenced  by  incomplete  conditional  equations,  k  least  ,  a 

squares  adjustment,  however,  trust  be  baaed  on  the  rigorous  conditional  equations 

designated  by  drolea  and  thin  lines  in  Figure  6* 


FIGURE  7 


Tha  corresponding  system  of  observational  equations  can  be  written  in  " 
analogy  to  formula  (2o)  as  follows i 

V  ..  ..  •  '  '  ..  :  •  :  ..  " 


where  A^  aim  the  coefficient  matrices  of  the  residual  veotara  Vj  correspond* 

log  to  tho  1  control  points  involved*  It  is  seen  from  Figure  6  and  Figure  7, 
that  six  different  types  of  A  A  submatrioes  exist,  depending  upon  what  type 

of  oontrol  point  ia  used  and  whether  the  point  is  oonmon  to  two  or  three 
photographs* 

The  A  submatrioes  of  the  types  II4  and  in^  are  obviously  identical 
with  the  corresponding  expressions  in  formulas  (2$), 

It  follows  from  formula  (98)  that  the  geometric  conditions  in  a  strip  can 
be  simulated  by  a  set  of  conditional  equations  explicit  in  terms  of  the 
residual  vectors  V •  Any  number  of  any  type  of  control  point  situated  at  any 
portion  of  the  atrip  can  be  handled  readily.  In  addition  ary  nunber  of 
Independently  determined  orientation  parameters  can  bo  enforced  in  the  solution 
at  ary  time. 

To  shew  the  procedure  of  a  atrip  triangulation  based  on  a  least  squares 
solution,  let  us  consider  a  strip  of  5  photographs  with  2/3  overlap  and  a  pass 
point  distribution  as  shown  in  Figure  8. 

Twelve  points  are  common  to  two  photographs  and  16  points  are  common  to 
three  photographs.  Consequently,  we  have  (12*2  ♦  18»3)*£-156  residuals,  which 
are  arranged  in  A  submatrioes  according  to  Figure  7. 
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The  number  of  observational  equations  is  given  byi 
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Thue  .there  ere  77  observational  equations  for  n  •  6  ■  5  •  6  »  30  unknown 
perimeters. 

With  the  symbols  used  in  Figure  7,  the  A  and  B  aetrLoes  of  the  observa¬ 
tional  equations  are  schematically  shown  in  Figure  9*  The  A  matrix  consists 
of  diagonally  arranged  A  4  sub matrices.  The  final  set  of  normal  equations 
can  b?<  formed  by  adding  sets  of  partial  normal  equations  as  outlined  on  pages 
22  and'’ 23.  Thus  the  analytical  treatment  of  a  least  squares  solution  for  a 
strip  triangulation  is  reduced  to  the  problem  of  inverting  a  symmetrical  r  x  r 
matrix,  the  coefficient  matrix  of  the  jc nknown  orientation  parameters.  Sines 
the  general  case  deals  with  six  unknowh  parameters  for  each  photograph  we  have 
for  a  strip  of  n  photographs,  r  ■  on.  II  The  fo«sibility°cf  the  proposed 
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approach  to  the  problem  of  aerial  triangulaticn  appears,  therefore,  to  be 
based  upon  the  possibility  of  inverting  largo  matrices.  This  problem  will  be 
sirrpllfied  because  the  individual  &  matrices  In  formula  (98)  include  only 
coefficients  corresponding  t >  the  unknown  parameters  combined  in  the  individual 
vectors  A.  Thus  the  r  x  r  natrix  of  the  coefficients  i3  only  partially  filled, 
namely,  in  the  neighborhood  of  the  diagonal  as  shown  in  Figure  10.  This  fact 
suggests  the  application  of  the  relaxation  technique  of  Qauaa-Seidel  or  a 
similar  approach.  The  determination  of  the  roots  of  the  normal  equations  can 
be  based  upon  an  iterative  process  which  takes  place  as  a  subroutine  within 
each  cycle  of  iteration  of  the  general  solution. 

If  the  pass  points  are  selected  in  such  a  way  that  at  leastf  5>  relative 
pontrol  points  are  common  to  each  section  of  triple  overlap,  those  models 
formed  by  the  combination  of  photographs  taken  at  every  other  camera  position 
are  included  in  the 'conputntiona.  The  favorable  geometry  of  the  thus  extended 
base  line  Is  otherwise  obtainable  only  by  convergent  photography.  If  we 
assume  2/3  overlap  between,,  the  consecutive  photographs  (see  figure  8),  the 
triangulation  strip  consists  of1  a  series  of  sections  wifch  triple  overlap. 
Consequently,  as  a  general  rule,  we  will  obtain  about  3  times  as  many  observa¬ 
tional  equations  as  we  have  relative -control" points.  In  addition  to  having  , 
favorable  geometry,  we  have  a  redundant  nunberof  points  which  are  used  to 
distribute  the  errors  of  observation  and  therefore  increase  the  accuracy  and 
re liability  of  the  result.  Photographs  taken  witty  lenses  of  hyper-wide  angle 
of  view  (120°),  which  are  expected  to  be  available'  in  the  future,  should  provide 
the  necessary  econony  in  a  strip  with  2/3  overlap. 

The  basic  idea  underlying  the  solution  presented  for  a  strip  triangulation 
can  obviously  be  applied  to  tho  problem  of  a  least  squares  solution  for  a 
block  adjustment.  Depend!  upon  the  degree  of  side  lap,  the  pattern  of  the 
A  subraatrioea  will  vary  according  to  the  change  in  the  number  of  photographs 
Which  have  in  common  the  images  of  certain  control  points.  However,  the 
matrix  of  the  final  system  of  observational  equations  will  again  consist  of  a 
sequence  of  diagonally  arranged  A  submatrices,  so  that  the  final  normal 
equation  system  may  be  formed  stepwise.  It  should  be  noted  that  the  applica¬ 
tion  of  formula  (20),  as  explained  on  page  h2 >  for  the  purpose  of  combining 
adjoining  models  can  now  be  used  advantageous!}'  because  in  this  case  the  rjpr^V,, 
to  bt^contoined  do  not  lie  in  or  close  to  a  plane.  ’  fZ  /:/4„ 


VI.  CONClDTi  NG  REMARKS 

The  feasibility  of  the  presented  analytical  treatment  of  a  photogrammetric 
evaluation  problem  depends  upon  the  availability  of  electronic  conputing 
macliines  with  large  computing  capacity  and  large  storage  facilities.  The  high 
speed  with  wliich  such  computers  are  able  to  perform  a  multitude  of  arithmetic 
operations  9ugge3t3  the  practicability  of  iterative  processes.  Therefore,  the 
least  square"  solution  is  advantageously  based  on  equations  which  yield  first 
order  differential  corrections  to  approximate  values  of  the  unknown  parameters,* 
the  reduction  bein^  repeated  until  the  solution  has  converged  to  a  pre-estabif$& 
ed  accuracy  level,  Tne  large  computing  capacity  rakes  it  feasible  to  design 


the  solution  for  the  most  general  case,  thus  eliminating  the  need  for  speotal 
solutions.  In  addition,  the  overall  accuracy  can  he  increased  by 
deteraining  the  solution  through  a  redundant  number  of  observations.  Any  type 
of  given  information  in  the  form  of  pass  points  or  orientation  elements  can 
be  rerdily  introduced,  ' 

Numerous  hand  oomputed  examples  have  been  carried  out  for  the  method 
described  in  this  report.  At  present  a  universal  code  is  being  prepared  at 
the  Ballistic  Rerearch  Laboratories  for  tha  problem  of  orienting  individual 
models  s*td  triai^ulatirg  the  spatial  coordinates  of  Individual  points  on  suoh 
models,  " 


It  is  assumed  that  suoh  Sjj  codo  together  with  an  automatically  recording 
stereoconparator  will  provide  a  solution  of  maxiniujn  preoioion’and  eoontay. 

It  appears  that  suoh  a  solution  is  applicable  to  all  photogrammetrlo  problems 
in  which  either  the  elements  of  orientation  or  the  spatial  coordinates  of 
numerous  target  points  must  be  determined.  .  a 
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